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Abstract 

The general method to obtain solutions of the Maxwellian equa- 
tions from scalar representatives is developed and applied to the diffrac- 
tion of electromagnetic waves. Kirchhoff 's integral is modified to pro- 
vide explicit expressions for these representatives. The respective in- 
tegrals are then evaluated using the method of stationary phase in two 
dimensions. Hitherto unknown formulae for the polarization appear 
as well as for imaging by diffraction. Ready-to-use formulae describing 
Fresnel diffraction behind a round stop are presented. 
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PACS numbers: 41.; 42.; 42.25.Fx; 42.25.Ja 

(The final version appeared in Z. Naturforsch. 65a, 1-24 (2010).) 

1 Sad state of theory 

The first impetus to this work came when your author read an altogether 
reasonable preprint on the Casimir effect. Its author, though obviously being 
a skilled physicist, wrestled heftily to establish the normal modes of the 
electromagnetic field between two metallic plates. The modified preprint is 
now published [lj. 

The second impetus came when your author lectured on theoretical op- 
tics. Diffraction of light belongs to the canon. He explained to the students 
what he had found in the treatise of Nobelist Max Born [2] and realized only 
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after a few weeks how bad it is. Born's book and its enhancement [3J are sub- 
titled "Electromagnetic Theory of Propagation, Interference and Diffraction 
of Light" . A truthful subtitle should rather be "A Wrong Theory of Sound" , 
at least with respect to diffraction. 

Nobody, neither physicist nor mathematician, seems to be able to system- 
atically solve partial differential equations for vector fields. But the Maxwell 
equations, foundation of optics and quantum electrodynamics, are just equa- 
tions of that kin. As a consequence, the theory of diffraction is still a grave- 
yard of intellectual atrocities. 

"Graveyard" in German means "Kirchhof" . It is Gustav Kirchhoff who 
is held responsible for the theory stolid till now [1],[5],[6]. Even in most 
recent textbooks Kirchhoff 's mistakes are recited as the theory of diffraction 
[I] , [E] , [H] , PH] , see [UJ for a historic review. 

Kirchhoff was aware of Maxwell's work, but did not esteem it too much. 
He preferred a product made in Germany, namely the Helmholtz equation 

V^(r P ) + A;V(i- p ) = (1) 

wherein r p is to denote the place of the probe, V p the nabla operator acting 
on r p and k the wave number, ip is the world-reputed symbol of parapsychol- 
ogy. Nobody knows what it means. Kirchhoff baptized it "Verriickung eines 
Aethertheilchens" translatable as "insanity of an etheral piece" 0, p. 9]. 

Today Maxwell's equations are not despised that much. It is generally 
accepted that light is an electromagnetic wave. Nevertheless scientists, even 
highly qualified ones, stick to Kirchhoff 's doctrines. They observe sagaciously 
that in Cartesian coordinates one may extract from Maxwell's equations one 
Helmholtz equation for each and every Cartesian component of the electric 
field E and the magnetic one B [3, Chap. 3]. The value of this observation is 
zero. The components of the electromagnetic vectors are coupled through the 
Maxwellian equations for curls and divergences. Solving the six Helmholtz 
equations for the components of E and B independently generates grossly 
wrong results. Even worse: In problems where Cartesian coordinates do not 
suit, as diffraction by spheres, there is, for the components of E and B, no 
Helmholtz equation at all \1'2\ . 

The next row of atrocities was digged when Kirchhoff, following Hermann 
Helmholtz, considered an integral which, seemingly, solves the Helmholtz 
equation [6, p. 82] 

= 1 , , f gxpg^-j) _ e*p(a|r p -r|U 

47r J Jf \ \r p — r| |r p — r| / 

wherein the surface F is to divide the entire three-dimensional space in an 
inner and an outer part; the inner part is where the probe r p resides, r points 
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to a point on F; it incorporates the variables of integration as components. 
df is the respective element of the surface. d n symbolizes a differentiation in 
the direction of the outer normal on F. The integral is interpreted in a way 
that primary waves ip(r) travel through outer space until they strike F; there 
they excite secondary spherical waves exp(ik\r p — r|)/|r p — r| interfering to 
produce the wanted i/j(r p ). 

Arnold Sommerfeld has criticized that simultaneous fixing of boundary 
values for the function and its derivative as required in Kirchhoff's integral 
((21) causes contradictions p3],[II],|15j. But scientists answered and still an- 
swer, they would take for ip(r) a plane wave and differentiate it consistently, 
without contradiction. These scientists do not understand that diffraction 
is possible only if part of the primary wave is screened. How should one 
treat these parts of Fl Kirchhoff and his followers demand simultaneously 
ip{r) =0 and d n ip(r) = on the screened part of F and call it a black screen 
[6j p. 40]. Yet the Helmholtz equation is elliptic of second degree. It has no 
real characteristics. Hence it follows from the Cauchy Kovalevskaya Theo- 
rem that the only solution compatible with a black screen is a global zero. 
Kirchhoff's integral produces no better than order-of-magnitude results. 

The third row of atrocites was graven when Kirchhoff wanted to evaluate 
the integral ()2]), but could not do it exactly. Especially the function in the 
exponents seemed to be invincible. Kirchhoff did then what all physicists do 
when they lack ideas. He substituted for the invincible function its crippled 
Taylor expansion. The evaluation with linear terms in r only was called 
Fraunhofer diffraction, whereas truncation after quadratic terms in r was 
said to yield Fresnel diffraction [61 p. 86]. 

One should assume that this kind of Fresnel diffraction includes Fraun- 
hofer diffraction as the special case in which quadratic terms are negligible. 
But this is not true. The formulae derived in this way describe fundamentally 
different pattern. 

All these problems will be solved in the present article. 

The general method to decouple Maxwell's equations will be developed in 
section [21 The electromagnetic field is represented by three scalar functions a, 
b and c which obey separated differential equations. One may solve the sepa- 
rated equations for these mathematical auxiliaries one by one and afterwards 
calculate the physical fields E and B from the auxiliaries by straightforward 
differentiation, c is the familiar scalar potential, a and b are common factors 
in the components of the simple and the double vector potential, respec- 
tively. This will be demonstrated for all electrodynamics in homogeneous 
and isotropic materials, in particular for those with finite conductivity. The 
General Representation Theorem in section [21 constitutes the first main result 
in this article. 
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In electrodynamics, diffraction is just a simple special case for which the 
representatives a and b suffice. 

In all problems with partial differential equations, also in the theory of 
diffraction, suitable boundary values must be fixed. Instead of the inconsis- 
tent black screen we will put up a screen of perfect conductivity. This entails 
boundary conditions for the electrical field E. When these are converted for 
a and b, it turns out that the one representative obeys a homogeneous bound- 
ary condition of the first kind (aka Dirichlet), whereas the other is subject 
to a homogenous boundary condition of the second kind (aka Neumann), see 
section [3j 

To complete the mathematical definition of diffraction, initial values must 
be given. In section H] we will account for them using a Laplace transform. 
Then the transformed representatives a and b fulfill separated Helmholtz 
equations. This explains why Kirchhoff 's formulas are not entirely wrong. 

For the representation theorem the way the representatives are found 
does not matter. Expansions in terms of partial waves are, for example, 
feasible, but the aim of this article is to mend Kirchhoff's advances. Since 
the representatives a and b fulfill Helmholtz equations, one may derive for 
them integrals similar to (j2J). As Sommerfeld's criticism must be attended, 
the spherical wave will be replaced with genuine Green functions for boundary 
problems of the first and the second kind in section [5j 

Hence it might seem that, for the true theory of electromagnetic diffrac- 
tion, expense is double as two integrals must be computed. Yet if for the 
primary wave a spherical wave is appointed that proceeds from a source at 
r q , both a and b are determined by the same function of r p and r q , the only 
difference being that probe and source are interchanged. The equations (157]) 
through ( 16 Up display the second main result in this article. 

The invincible function in the integral, mentioned above, can be rewritten 
to have an intuitive meaning: It is the difference of lengths when the points 
r q and r p are connected either directly or via a point r on the screen. These 
very terms occur in the familiar triangle inequation. It follows from the 
properties of the triangle function that simple geometry determines the basics 
of diffraction. For example, the border of shadow can be derived from the 
zero of the triangle function. The Criterion of Light will be established in 
section [6j 

Augustin Fresnel introduced, for the description of diffraction, certain in- 
tegrals. For a realistic theory of diffraction these integrals are too clumpsy. 
It is more convienient to use instead the error function known to all statis- 
ticians. Although one needs, for diffraction, the error function of a complex 
argument, it behaves in a similar way as that of a real argument. Aside 
for greater simplicity, we get from the error function more results, namely 
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diffraction in dissipative materials, which is of paramount importance for 
practical purposes. The salient properties of the complex error function will 
be described in section [7] 

There are several methods to evaluate the integral (158|) . e.g. for great 
distances of the probe from the screen. In this article we will apply the 
method of stationary phase to find the integral in the limit of short waves. It 
is this what Kirchhoff wanted when he expanded his functions up to second 
powers. We will avoid Taylor's expansion and introduce instead new variables 
that map the triangle function exactly. How this has to be done for two- 
and higher- dimensional integrals with finite limits, was hitherto unknown. 
It requires a novel deliberation. The Principle of Utter Exhaust will be 
introduced in section [HI It constitutes the third main result in this article. 

Utter exhaust (l§Uj) will be applied to the corrected Kirchhoff integral 
( 158]) in section [9j It yields the fourth main result in this article, the Uni- 
versal Formula of Diffraction (1 1 5 p : The diffracted wave equals the primary 
wave times the complementary error function which accounts through its ar- 
gument for the specific shape of the edge; the argument is found by a purely 
algebraic calculation. The formula is universal in so far as it holds for all 
single diffracting edges. The diffraction by screens with several edges can be 
derived from it by mere superposition. The universal formula also holds for 
probes arbitrarily close to the screen. Therefore it can describe the transition 
from Fresnel to Fraunhofer diffraction. It also holds for dissipative materi- 
als; the argument of the error function also accounts for damping. Yet the 
formula (I105p does not describe diffraction of very long waves. Moreover we 
will notice, under certain conditions, its weakness for great distances of the 
probe from the screen. 

In this article, no graph of calculated fields will be shown. Detailed com- 
parison to or prediction of measurements is here out of scope. Nevertheless 
some impressive applications will be outlined. The first is diffraction by a 
straight edge in section [TO]. Diffraction means breaking beams asunder. Nev- 
ertheless a metallic half screen creates via diffraction an image of the source, 
that is focusing of beams. Astonishing as this result might appear, the same 
formula contains as a limiting case Sommerfeld's stringent solution describ- 
ing the diffraction of a plain wave by a half plane, see section [TTJ Thus 
the asymptotic methods developed here have more power than Sommerfeld's 
stringent integration of Maxwell's equations. 

From the diffraction by a single straight edge it is just a small step to 
the diffraction by a slit, namely a superposition. Nevertheless the ensuing 
formula describes both Fresnel and Fraunhofer diffraction and the transition 
between these regimes, see section [121 

In section [13] the universal formula is applied to diffraction by a circular 
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aperture. It results in simple formulae describing Fresnel diffraction behind 
a round stop, a device that is used in almost all optical instruments. 

In the concluding section [HI the possibly novel results of this work will 
be listed and a program what is to be done next be given. The article ends 
with a surprise. 

This work is a hommage to the natural mind. All calculations were done 
using a pencil and sheets of paper. Aid from computers was rejected. 

2 The representatives of electrodynamics 

The purpose of the representation theorem is to replace Maxwell's vector 
equations for the magnetic and electric fields B and E with separated differ- 
ential equations for scalar representatives a, b and c. As soon as the latter 
equations are solved, successively, one may determine the vector fields from 
the scalar representatives by straightforward differentiation. To understand 
and to prove the general representation theorem of electrodynamics, three 
requisites are needed. First, the 

Lemma of Triple Curl. The vector differential equation 

V x V x V x va(r, t) = -Vx D t va(r, t) (3) 
can be reduced to the scalar differential equation 

VV r , f ) = A^, f ) + {^ * ZZl «> 

if and only if the supporting vector field v is preformed as 

v = v + v x r (5) 

with arbitrary constants vo and v\ . D t symbolizes an operator possibly in- 
cluding differentiations with respect to time t, but definitely no differentiation 
with respect to space r. /(•,£) denotes a free function. 

D t is, for example, efidf + fiad t , see equation ([T3"]) below, with constants e, 
\i and a denoting dielectric constant, magnetic permeability and conductivity, 
respectively. /(•,£) is a so-called gauge, a function which can be chosen 
according to convenience. For present purposes /(•,£) = suffices. 

The lemma was proven in |16j . It was published in [T7] and is available 
now in a textbook [18]. It is useful for uncoupling all vector equations which 

^dn this section and the ensuing two, the index p at the position of the probe r p will 
be omitted. 
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describe physical phenomena in a homogeneous and isotropic space, e.g. in 
the theory of elasticity and liquidity. 

Second, it is assumed that the electromagetic field propagates in an 
isotropic and at least piecewise homogeneous medium. The constitutive re- 
lations D = eE and H = B//x, which relate the force-exerting fields E and 
B with the source-caused fields D and H, as well as Ohm's law 

j = aE , (6) 

which relates the electric field E with the electric current density j, will be 
mounted in Maxwell's equations from the very start. 

Third, as we want to solve initial- and boundary- value problems in time t 
and three-dimensional space r for the electric field E(r, t) and the magnetic 
field B(r, t), it is necessary to discriminate between charges and currents 
that are enforced from outside and those that come about through the free 
play between inside fields. The charge density p (r) at the initial time t = 
t Q is determined by reasons outside the considered system. In conducting 
materials it decays exponentially. The remainder p(r, t) is zero at t = t and 
thus determined by the density of the current because of continuity 

d t (p(r,t) +p (r)exp (?-{t - *))) = -V(j e (r,t) + j(r,t) + j D (r,t)). (7) 

The current density, in its turn, is partly generated by external sources 
j e (r, t). In conducting materials, moreover, the electric field drives inner 
currents j(r, t) +j (r, t) according to Ohm's law (J6]). The latter is the current 
caused by the decay of p (r). 

General Representation Theorem of Electrodynamics. Solutions of 



the Maxwellian equations 

V x B(r, t) = epd t E{r, t) + paE{r, t) + /ij e (r, t) (8) 

VxE(r,t) = -^B(r,t) (9) 

VE(r, t) = ~(>(r, t) + Po (r) exp {^{t - t)) ) (10) 

VB(r,t) = (11) 



are provided by the solutions of the differential equations for the scalar rep- 
resentatives a, b and c 

V 2 a(r, t) = epd^a(r, t) + pad t a(r, t) — p J j e (r, r) exp ^— (r — t)^j dr (12) 
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V 2 6(r, t) = e/i3 ( 2 6(r, t) + i^ad t b(r, t) (13) 

V 2 c(r) = -~p (r) (14) 

if the magnetic and electric Gelds B and E are computed from a, b and c 
according to 

B(r, t) = V x v( j + St) a(r, £) - V x V x v6(r, t) (15) 
E(r, t) = —V x V x va(r, t) + Vx v<9 t 6(r, t) 

-Vc(r)exp(~(* -t)) -i^j e (r,r)exp^(r-t))dr . (16) 

Tie supporting vector Geld v must be parallel to the density of the enforced 
current 

j e (r,t) = v Je (r,t) . (17) 

For the general three-dimensional density, this amounts to choosing three 
linearly independent supporting fields according to (0), introducing three 
representatives a and solving three equations / fl2j) . 

The equations (fT2i) through (TT7I) constitute the first major item of this 
article. 

The proof proceeds in three steps since Maxwell's equations constitute a 
linear system. We compose the general solution from particular ones. Let us 
begin with the extraction of nonhomogeneities. 

First step: The ansatz 

E(r,t) = -Vc(r)exp^(t -t)) B(r,t) = (18) 

is to extract the nonhomogeneity of (flQj) with p (r)- The ansatz satisfies all 
Maxwellian equations except the third ( TTOj) . The third yields the Poisson 
equation (I14p for the scalar potential c(r). Please find ansatz (fT8l linearly 
enclosed in the general representation formulae (1151) and (j!6j) . 
The Ohmian current caused by this electric field 

j (r,t) = -aVc(r)exp^(t -t)) (19) 

balances in (J2J) the term with the charge density p (r). We discard them 
both to go on with a simplified equation of continuity and observe that it 
still facilitates the elimination of p(r, t). 



8 



Hence, for the remaining nonhomogeneities of Maxwell's equations, we do 
not miss anything when we differentiate the remainder of ffTUl) with respect 
to time 

V($E(r, t) + ^E(r, t) + p e {r, t)) = 0. (20) 

The equation of continuity ([7]) was used to get rid of the charge density 
p(r,t). Ohm's law ([BJ was applied to eliminate the current density j(r, t). 
Equation (f20|) . however, is guaranteed if the first Maxwellian equation flH]) is 
fulfilled. To see this, one just has to take its divergence. 

We can completely forget about the third and fourth Maxwellians fTTOl - 
[TT]) when we represent the magnetic field B by a vector potential A(r, t), i.e. 
B(r, t) = V x A(r, t). This shall be done henceforth. One should, however, 
keep in mind that mere introduction of a vector potential does not help much. 
The differential equations for the vector potential are coupled in a similiar 
way as Maxwell's equations for the electric and magnetic fields. We must 
construct special vector potentials to obtain uncoupled differential equations. 

Second step: The ansatz 

B(r,*) = V x vo(r,f) (21) 

introduces, as announced, a special vector potential, A(r, t) = va(r,t), viz. 
a predetermined vector field times a free scalar function. 

It shall be used to extract the nonhomogeneity j e (r, t). The first Maxwellian 
(jHJ) can be written as 

d t E(r, t) + -E(r, t) + -j e (r, «) = -VxVx va(r, t) . (22) 
e e e/j, 

Integrating this equation with respect to time yields a representation of the 
electric field 



E(r,t) = —V x V x v / a(r,r)exp (~(r-t))dT 

EfJ, ho \£ J 

J* Je (r, r) exp ^(r - t)) dr . (23) 

Most people prefer differentiations over integrations. Therefore we redefine 
the representative 

a (r,t) = ^"*a(r,r)exp (^(r-t^jdr . (24) 
This transforms (|21|) to 

B(r,t) = Vxv^ + 9 ( Ja(r,t) (25) 
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and ([23]) to 



E(r,t) = —V x V x va(r,t) - - /' j e (r, r) exp f-(r - i)W . (26) 
e\i e Jto V s J 

Please find these terms enclosed in the representation formulae (1T5]) and (TL6]) . 

The only Maxwellian equation (]9]) not yet fulfilled produces after insertion 
of ( 125]) and (1261) the nonhomogeneous telegraph equation (fT2]) . To see this, 
one has to apply the lemma of triple curl 02]) • This is possible only if the 
current density is parallel to the supporting field, i.e. if condition ([17]) is 
fulfilled. 

Now that we have taken care of the enforced current j e (r, t), we may 
assume for the rest of the proof that only the current caused by the inner 
electric field via Ohm's law (EJ) remains. This is not a nonhomogeneity. The 
first Maxwellian ([S]) can be simplified 

V x B(r, t) = //<rE(r, t) + e/id t E(r, t) . (27) 

Third step: The ansatz 

E(r,t) = V x v/3(r,t) (28) 
inserted into the second Maxwellian produces the representation 

B(r, i) = -VxVxvf 0(r, r)dr . (29) 

Again, for calculational convenience we redefine 

6(r,f)= fp{r,r)&r (30) 

Jt 

to obtain from ( 128]) and (129]) the still missing terms in the general repre- 
sentation formulae ( 1T5]) and ffTB"]) . The only Maxwellian which is yet not 
satisfied, the first one ( 127]) . causes, after application of the lemma of triple 
curl ([3]), the condition f ]T3"]) . again a telegraph equation, this time, however, 
a homogeneous one. Q.E.D. 

The general representation theorem copes with almost everything dis- 
cussed in most monographies on electrodynamics [19] : electrostatics, mag- 
netostatics, also electric discarge, skin effects in conductors, propagation of 
electromagnetic waves in space, guides and resonating cavities, optics, metal- 
lic optics too, radiation from antennae and all kinds of electromagnetic scat- 
tering, especially Mie scattering, which appears in thoses monographies as an 
extremely difficult case. For the theory of supraconductivity, Ohm's law (jHJ) 
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must be replaced with London's equation, but this even simplifies the deriva- 
tion of a slightly modified representation theorem. Only the electrodynamics 
of nonisotropic materials cannot be tackled in this way. 

Vector potential and double vector potential, dubbed Hertzian vector, are 
known since long. The furthest reaching representation of electromagnetic 
fields in terms of these potentials was probably found by Max von Laue [20J. 
For instance, the term in (I16p with the integral over the current density, 
which bewilders saplings, was given by Laue, but without consideration of 
conductivity. However, Laue did not disentangle Maxwell's equations. Many 
scientists subject their vector potentials to so-called Coulomb or Lorentz 
gauges. These gauges are related to invariances, but initial and boundary 
data break them. Thus in initial- and boundary-value problems, these vector 
potentials mislead. What one must use instead are adaptive scalars times 
fixed vectors. One can imagine suitable vector potentials as scalars riding on 
prepared vector fields like trains ride on rails. Inspiring in this direction was 
Peter Debye's simplified solution of Mie's problem using a special Hertzian 
vector [2TJ . A remarkably complete list of scalar functions that are useful for 
the disentanglement of Maxwell's equations was presented by Meixner and 
Schafke [22J . Yet this list is valid only for free propagation of harmonic elec- 
tromagnetic waves. The general principle of representation, i.e. the lemma 
of triple curl, was obviously unknown to all these scientists. Later on elec- 
tronic computers spread and absorbed interest. So this gap in mathematical 
physics was filled only in 1985 [16J. 

The brightest indication that scientists do not understand the general 
principle of representation is the lack of a reasonable theory of electromag- 
netic diffraction. There were attempts, for example |23j,[24j,[25j, (26], to ac- 
count for the vector fields, but they produced after longish explications only 
approximations - if at all. Also, partial wave expansions are not helpful since 
they converge well only for long waves [22] . By contrast, the theory that will 
be developed in the following sections is straightforward, yields exact equa- 
tions and is easily applied to practical problems. Kirchhoff or his scholars 
had done this if they had only known the approach. 

The situation was similar in hydrodynamics. With the same methods 
as explained here it was possible to derive, for the first time, turbulence in 
pipes from the Navier-Stokes equation [27]. In 1989 your author published 
the prediction that pipe turbulence consists of transients |28j. It was verified 
experimentally in 2006 [29]. 
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3 Boundary values on perfect conductors 



While the representatives a, b and c obey separated differential equations, 
they are usually tied together in the conditions on the boundary- and initial 
values of B and E. Consider, for simplicity, the diffraction of electromagnetic 
waves where explicit consideration of sources is not necessary. We do not 
need a scalar potential c. a and b must solve only homogeneous telegraph 
equations, cf. ( TT2|) and ( lT3]h 

Waves are diffracted when impeded by a screen. The only way to get on 
with boundary conditions rather than with conditions of transition is to have 
the screen made of perfectly conducting material. Then, because of Ohm's 
law for cr — > oo, the electric field E cannot maintain, on the screen, any 
tangent component. The waves do not penetrate. 

Let t denote any vector tangent to the screen 5* not to be confounded 
with the time t. The boundary conditions follow from the representation 

— (V x V x va(r, t))t + d t (V x v6(r, t))t = for r e S . (31) 

These are two equations because there are two linearly independent tangen- 
tial vectors t on a two-dimensional boundary. All the more it is surprising 
that the two unknowns a and b can satisfy the next four equations 

(VxVxvfl(r,i))t = fl and (V x v6(r, t))t = for reS . (32) 

This is possible if the supporting field v is parallel or perpendicular to the 
diffracting screen. 

To prove this, define a local Cartesian coordinate system 

F = G X X + Gylj + G Z Z t = G X t X + Gyty (33) 

such that it unit vectors and G y be parallel to the screen, whereas e 2 
pierce it normally. Just the components t x and t y of the tangential vector 
t are arbitrary though constant. Consequently only x and y components of 
curl and double curl need to be considered if the boundary conditions (J32]) 
are to be satified. 

v = g z is according to ([5]) an admissible choice for the supporting vector 
field. The equations (132]) become 

(V x G z b)t = G x td y b — G y td x b = 
(V x V x G z a)t = G x td x d z a + G y td y d z a = . (34) 

The differentiations with respect to x and y are inner ones since the tangent 
vector t of equation fT32l) is spanned by and G y . Both d x b = and d y b = 
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are satisfied on the screen if b = 0. Equally, both d x d z a = and d y d z a = 
are satisfied on the screen if d z a = 0. 

Generally the supporting vector field v (J3]) is not constant. However, 
if it is normal to the surface, we can construe its length as a factor of the 
representative and repeat the preceding calculation. Hence the following 
theorem: 

Theorem on Boundary Conditions. A supporting vector field v normal 
to the surface S of a perfect conductor induces homogeneous boundary con- 
ditions of the first kind for the representative b, whereas the representative 
a must fulfill homogeneous boundary conditions of the second kind 

b(r,t) = and <9 n |v|a(r, t) = for reS. (35) 

Oppositely, a supporting vector field tangent to the surface induces homoge- 
neous boundary conditions of the first kind for the representative a, whereas 
the representative b must fulfill homogeneous boundary conditions of the 
second kind 

a(r, t) = and d n b(r, t) = for r e S . (36) 

d n denotes differentiation along the normal on S. 

For convenience of reference, your author bundled the essentials of this 
section in one theorem. Its second part still has to be proven. With the local 
coordinate system introduced above, v = e^ is according to (JS} an admissible 
choice for the supporting vector field, too. The equations fj32j) are now 

(V x e x b)t = e y td z b = 
(V x V x e x a)t = e x t(dl — efid^ — fiad t )a + e y td x d y a = . (37) 

For the first component of the latter equation the homogenous telegraph 
equation f[T2j) was exerted. Again, the differentiations with respect to x and 
y are inner ones. The same is true for the differentiations with respect to 
t because boundary conditions must hold for all times. Hence d z b = and 
a = on the screen. When this is written without Cartesian coordinates, the 
second part of the above theorem emerges. A correction with the length of 
the supporting vector is not necessary here because v as defined in (J5]) does 
not vary in the direction of the normal if it is perpendicular to that normal. 
Q.E.D. 

Screens do not enclose radiation. Much of it spreads in open space. Thus 
the boundary conditions need to be completed, namely by retardation 

a or 6(r, t) ~ ( /a ° r ^M_^/vgg) + 0{(r ^ e -«t/2e for | rj ^ ^ (3g) 
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meaning that waves trail away in nirvana and never return. The functions 
f a and ft, may depend on the direction of the radiation, but they depend 
on the distance |r| only through the compound argument |r| — t/y/sjl. 0(a) 
is E. Landau's order symbol to appraise neglected terms on the right-hand 
side, see e.g. j30j Section 1.1]. Notice: a denotes here the conductivity of the 
propagating medium, for example air. 

Hence we may have boundary conditions that do not couple the represen- 
tatives. Yet to profit from the theorem we need supporting vector fields that 
are everywhere either normal or tangent to the screen. According to (j^J) this 
can be achieved for four types of screens: for plane ones, for parts of spheres, 
for parts of cones and for parts of cylinders, i.e. for cylinders with arbitrary 
cross sections. In this article we will be busy enough to cope with diffraction 
by plane screens and will use the theorem with the tangent supporting vector 
field ([36D- 



4 Initial values transformed 

We account for arbitrary initial values using a Fourier or rather a Laplace 
transform. All fields are proportional to exp(— iut). 

B(r, t) = B fc (r) e"*"* E(r, t) = E fe (r) e~ iut (39) 
a(r, t) = a fc (r) e** b(r, t) = b k (r) e** (40) 



The telegraph equations (11211131) become thus Helmholtz equations 

V 2 a fc (r) + k 2 a k {r) = V 2 6 fc (r) + k 2 b k {r) = . (41) 

The wave number k depends on the frequency u according to 

k 2 = E[iuj 2 + i/iau . (42) 

One of the two parameters, u or k, can be chosen as real. In the latter case, 
a negative imaginary part of u describes fading with time t — > oo. This is 
expressed through 

-^- ! s +0 <^ < 43 > 



and corresponds to equation (1551) . In the first positive imaginary part 

of k describes attenuation in space as |r| — > oo. Modulus and phase of k can 
be read from 



k = \l e 2 ^i 2 u A + jj 2 a 2 uj 2 exp ( - arctan — ) . (44) 

ecu/ 
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The point to be made is that the phase of k varies only between and 7r/4. 
This will matter in the discussions of section [7J 

Finally we must translate retardation (|38|) into the language of Fourier 
transforms. Both representatives and bk behave as leaving spherical waves 

a or 6 fc (r) = F a or , exp ^ |r|) + (\r\~ 2 ) for |r| oo (45) 

|r| 

at which the real parts of k and u are supposed to carry the same sign. The 
scattering amplitudes F a and F b may depend on the wave number k and the 
direction of the radiation, but not on |r|. Though Sommerfeld preferred to 
write this as d\ r \ak{v) = ikauiv) + 0(|r|~ 2 ) etc., your author persists in calling 
(J45]) Sommerfeld's radiation condition. 

5 Kirchhoff 's theory corrected 

Altogether we found that, for a plane screen, the representatives cifc(r) and 
6fc(r) can be computed separately. Both satisfy Helmholtz equations (HTj) 
and exhibit the same behavior in infinity (l4"5j) . Yet on the screen cifc(r) must 
solve a boundary- value problem of the first kind, whereas bk(r) is subject to 
boundary conditions of the second kind, see equations (1551) . □ 

Modifications of Kirchhoff 's integral ([2]) inaugurated by Sommerfeld are 
useful to solve both boundary-value problems. Instead of the spherical wave 
exp(zfc|r p — r|)/|r p — r| flexible Green functions G(r p , r) 

Hr p ) = ^JJ p (g(t p , r) d n ^(r) - ^(r) d n G(r p , r)) df (46) 

are introduced. Kirchhoff 's formula remains valid if these Green functions 
fulfill the Helmholtz equation (Tj[|), if they have the same singularity as the 
spherical wave 

<9| rp _ r |G(r p , r) = -(r„ - r)" 2 + 0(\r p - r^ 1 ) for r p -> r (47) 

and respect the radiation condition as in equation (j4"5l) . 

In a boundary-value problem of the first kind, the values of the function 
■0(r) itself are known on the plane, but no information on the values of the 
derivative d n ip{y) is available beforehand. Hence we need a Green function 
that is zero on the plane lest the unknown values matter. We find from ( I46p 

= ^J j F d n G x {T p ,v) ^(r)d/ if G x (v p ,r) = for r e F . (48) 

2 From now on, the index p at the position of the probe r p is indispensable again as 
there are other locations which must be discriminated, viz. the position of the source r q 
and arbitrary points on the screen r. 
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On the other hand, if the values of the derivative are known, we must do 
without the function itself. For the solution of the boundary-value problem 
of the second kind we need a second Green function 

i/>(r p ) = ^JJ p G 2 (r P , r) 9^(r) df if d n G 2 (r p , r) = for r e F . (49) 

For F being a plane characterized in Cartesian coordinates (I33p by, say, z=0, 
both Green functions are easily found as a spherical wave around the point 
of measurement r p plus its mirrored image on the other side of the plane 

Gi or 2 (r p , r) = _ or + exp(»fc|r - r|) (g()) 

with 

Fp G x Xp -\- GyUp ~\- G z Zp F m G^Xp -I - GyUp Gz^p • (^^) 

Both functions own the correct singularity as in (14"?]) and respect the radiation 
condition (T45|) . 

Furthermore it shall be assumed that the primary wave is created by a 
point-like source at r q behind the plane 

^{v) = eXp f^- r ^ . (52) 

I 1 " r q\ 

This kind of primary wave is to be preferred over the usual plane wave 
because it respects a radiation condition of type ( 1451) and ensures thus the 
validity of the generalized Kirchhoff formula (146 p . Also it is closer to exper- 
iments as it is simpler to produce an approximate spherical wave than an 
approximate plain wave. The beloved plain wave can be obtained from the 
spherical wave by a straightforward limiting transition \r q \ — > oo. Radiation 
from higher multipoles can be derived by differentiation with respect to 
and the general primary wave by superposition. 

Writing down the integral for the boundary conditions of the first kind 
(148]) in Cartesian coordinates ends up with 



- ^ (53) 

and the integral for the boundary conditions of the second kind ( I49p becomes 
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with 



R p (x,y) = l r p - r(x,y) \ = J(x - x p ) 2 + (y - y p ) 2 + z\ 



R q {x, y) = \t{x, y)-r q \ = y/(x- x q ) 2 + (y - y q ) 2 + z\ (55) 



R pq = \r p - r q \ = yj (xp - x q ) 2 + (y p - y q ) 2 + (z p - z q ) 2 

and the triangle function 

A(x,y) = R p (x,y) + R q (x,y) - R pq (56) 

which measures the difference of lengths when the points r p and r q are either 
directy connected or via an arbitrary point r = e x x + e y y on the intermediate 
plane. 

It is not by accident that the ip's in the equations (1531) and (1511) are 
identical. The factors behind the spherical wave exp(ikR pq ) / R pq including 
the double integrals extending from — oo to +oo have both the value 1; a 
proof of this fact will be given in section [TTJ The formulae fl53l) and fl54|) 
are just different mathematical realizations of Huygens' principle: A wave 
propagating from a source point that strikes a plane excites there secondary 
waves interfering to reproduce the original wave. 

Diffraction happens only if a screen covers parts of the plane. Using the 
boundary conditions in the theorem with the tangent supporting field (136]) 
yields explicit expressions for the representatives 

/ \ exp(ikR pq ) . . exp(ikRp q ) 
Wp) = p JJ ( r P' r i) b k( r p) = 5 - n ( r n r p) ( 57 ) 

lt pq £i> P q 

with the double integral 

n(r r ) - -\ z p\ R pi [[ ik-l/Rp(S,y) d(x,y) jkA&r,),,, ( . g) 
2n JJoR 2 p (Z, V )R q (Z, V ) d(£, V ) e ^ dT] - (58) 

To be general enough for all applications, arbitrary variables £ and rj shall 
substitute the Cartesian ones x(£, rj) and ?/(£, rf). The functional determinant 
d(x, y) j d(£, rj) enters to transform the element of the surface df = dxdy. The 
integral // extends only over the aperture or diaphragm D in the surface F. 
The arbitrariness of £ and r\ shall be used for a simple description of the 
aperture in the way that freely varying £ with a fixed rj depicts an edge. The 
aperture is described, for example, by — oo < £ < oo and — r\_ < r\ < rj + . 

The enthusing result is that, for the full electromagnetic theory of diffrac- 
tion with all possible polarizations, one needs to evaluate only one integral. 
There is scarcely more work to be done than in the theory for scalar waves. 
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The addititional work consists in some elementary differentiations as pre- 
scribed in the representation formulae (I15H16P simplified for Fourier trans- 
forms as 

B fc (r p ) = f j - »w) V p x ta k (r p ) - V p x V p x tb k (r p ) (59) 
E fe (r p ) = — V p x V p x ta fc (r p ) - zwV p x t& fe (r p ) (60) 

with a supporting vector field v=t that consists of a constant tangent as in 
AMD- 

The equations floT|) through flBTJj) constitute the second major item of 
this article. The integral ( |58|) can be evaluated using various techniques. 
Your author will discuss in the following sections only one: the method of 
stationary phase applicable for short waves with §lk — > oo and > as 
conditions on the real and imaginary parts of the wave number k, respectively. 

While equations ( 15"TH58"|) solve mathematically well-posed boundary-value 
problems without any error, they do not describe physical diffraction exactly. 
When a conducting screen diffracts a wave, reflection cannot be avoided. 
Most of the reflected wave stays in outer space, z p < 0, but it is also diffracted. 
A tiny fraction of reflection invades inner space z p > 0. We will learn to 
handle this in section [TT1 



6 Light and shadow 

The triangle function rules the diffraction of short waves. In its definition f[5T)j) 
interest was focused on the point of the screen r. Yet the triangle function 
also depends on the points of probe r p and source r q 

A(r, r p , r q ) = \r p - r| + |r - r q \ - \r p - r q \ . (61) 

Your author decided to locate 

screen r = e x x + e y y + e z z at z = , 

probe r p = e x x p + e y y p + e z z p at z p > , (62) 

source r q = e x x q + e y y q + e z z q at z q < . 

The triangle function is positive except at that point r = r s = e x x s + e y y s 
where the straight line between probe and source pierces the plane 

Xs + iys= ^ + ^)-^v + m) (63) 
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Complex notation is preferred because it eases transformation to other co- 
ordinate systems, see below. The point r s is the location of the absolute 
minimum of the triangle function, and all the more it is a stationary point. 

The consideration holds for fixed points of source and probe. However, 
if only r q is fixed, whereas r p varies while r = r s slides on the edge of the 
screen, then 

A(r, r p , r q ) = if r on the edge (64) 

determines as function of r p a surface, namely the border of shadow. 

For calculating diffraction, we will need the root of the triangle function. 
Your author utilizes the ambiguity of the root to demand 



+| JA(r,r p ,r g )| if r p in the shadow, 

A(r,r p ,r q ) = ^ ^ (65) 

A(r, r p , r q )\ if r p in the light. 



It is cogent to assign different signs to the dark and the bright if the function 
is to be different iable. The triangle function is an analytic function which 
depends, around its minimum defined by (16~4"1) . quadratically on its variables. 
Therefore omitting the signs in (165]) would induce a similar discontinuity as 
in the assignment = \x\. The absolute assignment of the sign, on the 
contrary, is arbitrary since diffraction either by a screen or its complement 
is equal; remember Babinet's principle [2J § 47]. 

To decide where there is light or shadow, a handy criterion is needed. 
There is light on the probe if the screen does not impede the straight con- 
nection between source r q and probe r p . Thus the positive sign in f )65|) is to 
be taken if the stationary point x s + iy s in f )63|) misses the aperture. 



Criterion of Light. When one uses transformed coordinates £,77 adapted 
to the screen such that the aperture is described by r\_ < n < r] + while £ 
varies freely, the negative sign of triangle function 's root lfU5\) has to be taken 
if 

i]- < rj s < r] + . (66) 
r) s is calculated from the stationary point / f63j) via coordinate transformation. 



For example on a screen with a circular aperture, cylindrical coordinates 
p, <f, z — suit. The aperture is defined by p < p with p as the radius of 
the stop, while ip varies freely. The transformation to cylindrical coordinates 
is facilitated through 

x s + iy s = p s e %ipa x p + iy p = p p e llfiv x q + iy q = p q e %Lpq (67) 
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with the p's as axial distances and the <p's as azimuthal angles. Insertion into 
dHS]) produces two equivalent formulas for the axial distance 



Ps 



yj {Z p pq ~ Zqp p ) 2 + 4z p Zqp p p q S 2 



\J (ZpPq + Z q p p ) 2 - Az p Z q PpP q C 2 



(68) 



with the abbreviations 



C = cos 



2 



S = sin 



<Pp ~ ^Pq 

2 



(69) 



p s > po is thus the criterion for the domain of shadow, i.e. for the positive 
sign in ( )65|) . 

Instead of applying elementary geometry, as was done in this section, 
one may calculate the stationary point £ s , r] s by simultaneous solution of the 
equations 



%A(£,77)=0 d v A(Z, v ) = Z = £ s r) = r] s . (70) 



The result is, of course, the same as that given in fl63j) with subsequent trans- 
formation of coordinates, but the computational effort is larger. Your author 
displays the equations flTOj) only to ease comprehension of the astonishing 
equation ( 155]) which will appear in section [5J 

7 Using the error function for diffraction 

In the theory of diffraction, Fresnel integrals 



are still custom, but the error function or rather the complementary error 
function is handier [3Tj . [32] 



All statisticians become perfect opticians when they are willing to handle 
their favorite function with complex argument. 




(71) 




(72) 
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The error function comprises the Fresnel integrals in a similar fashion as 
the exponential function contains sine and cosine 

eric(V^iz) = 1 - V^iV2 {{C x {z) + iS x {z)) . (73) 

If z is assumed as real, \f—iz varies on the second main diagonal of the 
complex plane since \f—% = (1 — i)/y/2 = exp(— ire /A). The features known 
for real argument remain if the complex argument is enclosed between the 
first and the second main diagonals of the complex plane; expressed by a 
relation between imaginary and real parts: \Qz\ < \$tz\. This is exactly what 
we need for optics, see (jHJ). For large negative real parts Viz — > — oo of the 
argument the complementary error function starts at the value 2, assumes 
the value 1 at the origin z — 0, and attains the value for large positive 
real parts Vtz — > +oo. In the crudest approximation, one may think of the 
complementary error function as 2 for negative arguments and for positive 
ones. It is a switch. 

The asymptotic expansion of the error function familiar on the real axis 
remains valid in the wedge between the main diagonals which contains the 
real axis 

ertctz) = —=r- e~ z \l + 0(\z\- 2 ) for Viz -> +oo . (74) 

y/TTZ 

The asymptotic expansion on the other side of the complex plane Viz < 
follows from 

erfc(z) = 2 — erfc(— z) . (75) 

The only refinement due to complexity is that the complementary error 
function decreases monotonously on the real axis, whereas it takes complex 
values and both real and imaginary parts oscillate when the argument be- 
comes complex. 

8 The method of stationary phase for two- 
dimensional integrals 

Let us recall the asymptotic calculation of one-dimensional integrals 

I(k,rj^,T] + ) = r + A{ V ) e ikAiv) dr] (76) 

for $lk — > oo with ^sk > 0. We assume that the real function A (77) > is 
stationary for some 77 = r] s , i.e. 

Mv) = ^ L (v-Vs) 2 + 0((r ] -r ]s f). (77) 
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The indices at A indicate that the function be differentiated twice and the 
result evaluated at 77 = rj s . 

The familiar approach in the method of stationary phase is to introduce 
a new variable 

v = ^(v~Vs) (78) 

and to forget the higher-order terms 0((i] — rj s ) 3 ) in equation (1771) . The 
function 



k 



5(v) ~ \l— e tkv for m ->• 00, %k > (79) 

may be construed as a representation of Dirac's delta function. Thus the 
integral ( 176]) yields 



/ 2m 

I(k, -00, +00) ~ J— A(r) s ) . (80) 

The amplitude A{rf) appears as a constant. 

If the limits of integration r]± are ±00, this is correct, but for finite limits 
the local approximation (ITHj) induces systematic errors. What we have to use 
instead is the global map 77 — » v 



To rewrite the integral from the variable 77 to the variable v, we must calculate 
the differential drj = (di]/dv)dv = (dv/dri)~ x dv. Because of (179|) it is sufficient 
to know the value of the differential for v = corresponding to 77 = rj s . Thus 
the value of the differential of the global map (IHTj) to be used in the integral is 
the same as the differential of the local approximation (178 p . The peculiarity 
of the global map appears only in the limits: 



/ 2-ki erfc J —ikA (r?_ ) - erfc J —ikA (rj + ) 

i(k,v-,v + ) - A( Vs )—« \LL- Y \J±i (82) 

While the preceding is not familiar, it is known [301 Section 2.9]. Yet in 
the theory of diffraction one needs to evaluate two-dimensional integrals 

/■77-j- poo 

II(k, 77_, n + )= / A({, rf) e lhA ^) dtd V . (83) 

J T)— J— OO 

Again it is assumed that the real function A(£, 77) > is stationary at £ s ,r/ s 

= ^ (£ - is? + A s?3 - 6)(r/ - %) + ^ (77 - T7,) 2 

+0(K-6| 3 + (|r7-r/ s | 3 ) . (84) 
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It seems to be a hitherto unsolved problem to find a suitable two-dimensional 
global map £, 77 — >■ u, v such that 

u 2 + v 2 = A(£,r/) . (85) 

Here is the solution: The map is 



u = JA(Z, V )-A S ( V ) (86) 



v = y/A.(ji). (87) 
The function A s (r/) is determined by the 

Principle of Utter Exhaust. Eliminate £ from tie derivative 

d^, v )=0 £ = S.fa) (88) 

to find the exhausting dependence £ = S s (r/). Insert the exhausting depen- 
dence into the function A (£,77) to obtain the exhausting function 

A s ( V )^A(E s ( V ),v) ■ (89) 
The integral j83|) is, for — >■ 00 and 3fc > 0, asymptotically equal to 



. 2yri A(£ s ,r/ S ) eiicJ -zkA 8 (ri-.) - erfcJ -ikA 8 (rj + ) 
II(k,r]-,r) + ) ~ —7== == 5 



(90) 



Utter exhaust follows from the following indispensable requirements: 

£ = £ 8 , 77 = ?7 S be mapped to it = 0, u = , (91) 
it, v be real for all £, 77 , (92) 
f 2 = /(//) be a function of 77 only . (93) 

We need the third requirement (|93|) to retain the simplicity of the limits in 
the integral II fl53"j) when mapping £, 77 to w, v. Because of equation flSBI) and 
the requirement (192]) the function f(rf) must never exceed A(£, 77) whatever 
value £ takes. Nevertheless for every 77 there must exist £ = ^(77) such that 



equality is reached: f(rj) = A(H s (t7), 77). Otherwise tj = y A(£, 77) — 7(77) 
cannot take 0, a contradiction to the requirement (191 p . In other words, /(r?) 
and A(£, 77) coincide at that value of £ = S s (?7) where A(£, 77) gets stationary 
with respect to £. Hence S s (t7) is determined by elimination of £ in the 
condition 
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The final formula ( 19"U|) can be understood when compared to formula 



The double integration generates the factor wi/k, the square of ^iri/k. 

The determinant in 2/ A^A m — A|^ of the map £, rj — > u,v replaces the 

differential in ^2/ A m of the map rj — > v. Q.E.D. 

The principle was named as of utter exhaust because the function A s (rj) 
is the largest possible function of 77 only that fulfills A s (rj) < A(£, 77); it takes 
from the two- variable function as much as a one- variable function can afford, 
cf. equation flSol) . 

A further intuitive interpretation follows from a comparison of equation 
(I88p with the equations (ITUl) . The stationary point is where the function 
A(£, 77) takes its absolute extremum. On the exhausting dependence, by con- 
trast, A (£,77) is extremized only with respect to the one variable £, whereas 
the other variable 77 is fixed. In optics, A (£,77) essentially measures the 
distance between points. The absolutely shortest connection is of course a 
straight line. But the exhausting function A s (t7) measures a conditionally 
shortest distance, namely if the connecting line is forced to touch the edge 
of the diffracting screen. One can determine the exhausting function experi- 
mentally using a ribbon of rubber and a lubricated model of the edge. 

The difficult part of utter exhaust is the elimination according to equation 
( }88]h It is therefore gratifying to possess linear approximations of equations 
(I86H87P . similar to the one-dimensional case (1781) . These approximations are 



//* (£-&) + -^=(17 - 17.) (94) 



\ 2A e{ 

They were found by utter exhaust (|88p applied to the quadratic terms on the 
right-hand side of (1841) : searching principal axes of the ellipse is not a good 
idea. For the integral 77 (1551) , equation (1901) can be used when the exhausting 
function is approximated as 

A s (r/ +0 r-) - A ^ A ;;" A|t? (r7 + or- - Vsf . (96) 
2A ?? 

The relation to the border of shadow discussed in section [6] appears here at 
first sight. 

The method of stationary phase is sometimes blamed as not being math- 
ematical stringent. It is argued that certain integrals do not converge and 
thus certain errors cannot be estimated. The criticism does not apply here. 
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Guided by mother nature, we made a theory where the parameter k of asymp- 
toticity has a positive imaginary part $sk > 0. This guarantees the conver- 
gence of those integrals. We can even write the complete asymptotic series. 

Augustin Fresnel invented the zone construction in 1816 to prove that 
light travels as a wave, but this was just a semi-quantitive idea [33j P-247], 
[3J § 8.2]. Its mathematical solution is the principle of utter exhaust presented 
only now. 



9 The universal formula of diffraction 

To make use of stationary phase for diffraction, the amplitude in the integral 

am 



2tt Rl(i s ,T] s )R q (i s ,r] s ) d(£,7i) 
is to be evaluated at the point of stationarity £, s ,Vs an d multiplied by the 



factor 27ri/kJ A^A m — A^ of equation fl90|) . 

Trivially, R p and R q are just the two fractions of R pq 

Rp{is,Vs) = ^ Rpq Rq(£s,Vs) — ^—Rpq- (98) 

Zp Zq Zp Zq 

Distances as they are, they do not depend on the coordinate system. This 
is different for the determinant A^A m — A|„. Let us calculate it first in 
Cartesian coordinates. Differentiating the triangle function ( 156]) twice and 
evaluating it at the stationary point gives 

A _ ( Z P ~ Zq) 2 ((y P ~ Vqf + ( Z P - Z q ) 2 ) 

—ZpZqltpq 

^ _ ( Z P ~ Z q) 2 (( X p ~ X q) 2 + ( Z P ~ Z q) 2 ) ^gg^ 
— ZpZqRpq 

^ _ ~( z p ~ z q) ( x p ~ x q)(yp ~ Vq) 

%n R ,,,, 



Consequently 



"p^q^'-pq 



, 2 _ (%> ^<?) 6 



A A - A = — — (100) 

^xx^yy ^ xy 2 7 2R4 " 

^p^q h pq 

The root of this determinant is transformed multiplying it by the functional 
determinant taken also at the point of stationarity 

d(x,y) 



\t\v - Aj v = ^A xx A yy - A% . (101) 
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Additional terms do not occur. They would contain first derivatives of the 
triangle function which are zero at the stationary point, see conditions (1701) . 
The result is astonishing: 

2ixi A(£ a ,7) 3 ) _ 1 , 



kjA^A m - Aj v 



Oik' 1 ) (102) 



though it should be noted that it is -1 if the functional determinant is nega- 
tive, i.e. if the coordinates £,77 form a left-handed system. The final formulae 
for the representatives derived from (190|) are thus simple and identical 



atfo) = expg^) erfcy-ttA.t,-) -ericy-flfcA.^) +Q(t . 1) ^ 

Rpq 2 



exv(ikR„„) erfc-i/— ikA s (ri) — erfc-i/— ikA s (ri + ) 
&fe ( rp ) = Y ^WJ V tllli + OCA;- 1 ) . (104) 

Rpq 2 

Further simplifaction is possible. If the upper edge is removed, i.e. r/ + — > 
+00, the exhausting function is supposed to increase beyond measure A s (r/ + ) 
+00 and, according to the asymptotics (174"j) . the formulae (jl03til04j) are re- 
duced to the 

Universal Formula of Diffraction. The diffracted wave is the primary 
wave times a universal function describing the change from light to shadow. 



exp(ikR pq ) erfty-ifcAs nfu -u / inK N 
a or b k {r p ) = h 0(k ) . (105) 

The geometry of the diffracting edge enters just through the argument of 
that function. This argument is calculated from the exhausting function 
( 189)) evaluated at the position of the edge: A s = A s (r/_). 

From here one may return to the apparently more general formulae above 
due to the linearity of Maxwell's equations and a superposition of their so- 
lutions. So A s may be understood as an abbreviation for A s (r/_) or A s (r/ + ) 
as required. By the way linearity: Factors on the right-hand sides that do 
not depend on r p are always free. Your author omitted them. This is all the 
more tolerable since the representatives get effective only if taken together 
with the supporting vector field t in equations (1591) and (1601) . According to 
(J33]) the tangent vector t contains the components t x and t y . They can be 
chosen to adjust strength and polarization of the primary wave according to 
experimental givens. 

Formula fl 1 5 [) constitutes the main achievement of this article. 
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One may check the accuracy of (11051) by insertion into the Helmholtz 
equations (jHJ). One finds that (11051) satisfies these equations in the order of 
k 2 identically, but the terms proportional to k 3 ^ 2 cancel only if the exhausting 
function A s satisfies an equation of eikonal type. 

(V P A S ) 2 = -2^— ^V P A S . (106) 

Indeed, the fulfillment of this equation follows from the definition of the 
triangle function f lST]) and the condition of utter exhaust fl88l) . Hence the 
largest erroneous term is O(k). But since the order of the Helmholtz operator 
is k 2 , the relative error is 0(k~ l ), as expected. 

Checking the fulfillment of the boundary conditions yields at first sight a 
worse result, namely a relative error of 0(k~ l l 2 ). We see this from the asymp- 
totic expansion of the complementary error function (TM|) . This asymptotic 
expression applies because, directly behind the screen, the probe resides in 
the shadow where there is, according to (I6"5"|) . ^/~A~ S > 0. Thus the argument 
of the error function y/—ikA s lies in the right-hand side of the complex plane. 
Hence 

= exp^ exp^ + 
K pg 2,\J— iixkA s 



instead of an exact zero as required in the boundary conditions fl36j) . A 
similar result holds for the representative bk(r p ). To verify its boundary 
condition, one must calculate its normal derivative, i.e. its derivative with 
respect to z. While this derivative generally is O(k), its value on the screen 
is 0{k 1 ^ 2 ) such that the relative error is 0(k~ l l 2 ). 

At the second view, both results are better than expected. First, in 
Kirchhoff's scalar theory of diffraction, the probe must not approach the 
screen since Kirchhoff's functions become singular [6j p. 9]. In the present 
theory, all expressions stay regular. This facilitates, as we will in section 
[TTj a simple trick to suppress errors on the boundaries. Second, formula 
(I105p . even as it is now, satisfies the boundary conditions not exactly, but 
accurately. The reason is the product kA s in the denominator of (I107p . The 
exhausting function A s measures, in a slightly unusual way, the distance of 
the probe from the border of shadow. The distance between that border and 
the screen is the largest possible distance available in a diffracting system. 
Hence for short waves, both the wave number k and the distance A s are big 
making (I107p negligibly small. 

For applications it is worthwile to calculate from fl 1 5 [) the field strengths 
according to (159]) and (160]) . For simplicity we let dfc(r p ) = 0. Then the 
electric field Ejt(r p ) stems only from the simple curl in (I6T)]) and the magnetic 
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field B fc (r p ) only from the double curl in (J5Uj) . 



R 



(r p — r q ) x t erfci/— ^A s 



exp(ikR pq ) 



R 



cok + 



exp(zfcA a W — 



+ 0(fc) 



:io8i 



Tj / \ (r p - r,)((r p - r,)t) - tR 2 ericV~ikA s 2 
Jjfc^rpj — K + 

exp(ikR pq ) R l pq 2 

V p A s ((r p -r 9 )t) + (r p -r g + J R w V p A s )(V p A s t) /ifc3 
— j= exp{tkA s ) \ - — 

+0(k) . (109) 

The ubiquitous primary wave was drawn to the left-hand sides to direct 
attention to the nontrivial terms on the right-hand sides. 

We learn from these equations that diffraction of electromagnetic waves 
without polarization does not exist. Yet the contributions to the right-hand 
sides with the complementary error function describe just the vector proper- 
ties of the primary wave. The changes of polarization effected by the screen 
are described by the contributions with the exponential function. The con- 
tributions with the error function are proportional to k 2 since u = 0(k), 
cf. equation fH3|) . whereas those with the exponential function are propor- 
tional only to A; 3 / 2 . Nevertheless one should not be deluded that polarization 
by diffraction is always a lower-order effect. In the shadow the error func- 
tion decreases dramatically. It ceases to be 0(1) and weakens to be only 
0(k~ 1 ^ 2 ) as indicated in equation (11071) . In the shadow both contributions 
to the right-hand sides of fllQ8j) and (I109p reach comparable sizes. We must 
expect hitherto unseen effects. 

All terms in (llOSp and (11091) remain finite in the entire domain of solution 
except, of course, immediately at the edge where Bremmer's effect takes 
place. The singularities on the border of shadow caused by \/Aj in the 
denominators are cancelled by the derivatives V P A S in the nominators since 
A s varies quadratically in the vicinity of the border. 

For radiation with low frequencies, e.g. microwaves, the electromagnetic 
fields can be measured directly. However, the Maxwell equations are real 
equations for real observables. Thus we must extract from the fields (llOSp - 
(11091) their real parts prior to comparison with experimental data. In optics 
with its much higher frequencies, the classic measurements are all calorimet- 
ric ones. One collects the energy flux impinging on a given surface for times 
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much longer than inverses of these frequencies. The observable is here the 
time average of the Pointing vector S 



1 r dt 1 

S(r p ) = lim - / &E(r p ,t) x &B(r p ,t)- = — K(E fc (r p ) x B*(r p )) . (110) 

The relation between time-dependent and time-independent fields is as 
defined in equations ( 159"|) . The asterisk denotes complex conjugation; one 
can affix it either to the magnetic or to the electric field without affecting 
the result. The a priori decomposition of the complex fields (110811109]) is not 
needed here. 



10 Diffraction by a straight edge 

The first application is, of course, diffraction by a straight edge. For its 
description, the Cartesian coordinate system is best. Hence £ = x, r] — y 
and 

— oo < x < +oo — oo < y < y_ = const z = (HI) 

defines the screen. y_ is the location of the edge. 

The exhausting dependence is found from equations (|88|) and (1561) 



d x Rp(x,y) = -d x R q (x,y) . (112) 

Squaring this equation and cancelling identical terms on both sides results 
in a quadratic equation for x 

(x p - x) 2 ((y q - y) 2 + z 2 q )) = (x q - x) 2 ((y p - yf + z 2 p )) . (113) 

y must be considered as given. One solution of this equation is a mathemat- 
ical ghost. It appears because the original equation (11121) was squared to get 
rid of the roots. The other solution, however, does solve the original. It is the 
searched-for exhausting dependence £ = 5 s (t;) fl88|) in Cartesian coordinates: 

x P J (y q - y) 2 + z 2 + xj (y p - y) 2 + z 2 p 

x = X s {y) = — -, . — . (114) 

V (y q ~ y) 2 + z\ + yj (y p - y) 2 + z 2 



Inserting this dependence in the triangle function (I56p gives the exhausting 
function of the straight edge 

A.(y) = \J{x P ~ x q ) 2 + (yf(y p - y) 2 + z 2 + y/(y q - y) 2 + z 2 f 
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We take this function at the edge y = y~ and need to determine the sign 



of a/AJ = yA s i.e. where there is shadow or light. Shadow prevails, 
according to the criterion (155]) . when y s < y_. With y s from equation (153"]) 
we find 

y P <V- + (y--y q )^ E - . (116) 



Hence equation ( 155]) reads here 



= J +lv / Al| if y P <y-- (y- - y q )z P /z q with z p > o, , n7 * 

* s \ — \\fK~ s \ elsewhere. 

What remains to be done is to insert this in the universal formula of diffrac- 
tion ( IT05]) . 



11 Imaging by diffraction 

We have now enough experience to perform the simple trick that was an- 
nounced in section |9] To improve the fulfillment of boundary conditions, one 
inserts a mirrored picture of the source in the same way as Sommerfeld did 
when he calculated suitable Green functions, see section If r q is the posi- 
tion of the source, r m has the same coordinates x q and y q , but the opposite 
value of z q as explained in ( 15T]) . Let the distances R pq and R pm be defined 
according to ( 155]) : 

Rpg = \l (xp ~ x q ) 2 + (y P - y q ) 2 + (z P - z q ) 2 
Rpm = yj (x p - x q ) 2 + (y p - y q ) 2 + (z p + z q ) 2 . (118) 

Hence original and mirrored distance become equal, \im.R pq = limR pm , if 
the probe approaches the screen, z p — > 0; it does not matter if the approach 
takes place in inner or outer space, z p > or z p < 0. 

We must, however, be careful when we introduce the mirrored exhausting 
function. In the definition of the original exhausting function taken at the 
edge y = y_ 

A sg = yj{x p - x q ) 2 + {^y~^ y.) 2 + z 2 + ^y~ q -y-) 2 + z 2 f 

-yj(x p - Xq) 2 + (y p - y q ) 2 + (z p - z q ) 2 

+\\f^q\ if y P <y-~ (y- - y q )zp/z q with z p > o, 

I V ^/ 

— \JA S q\ elsewhere. 
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the assignment of signs of the root belongs to the definition, cf. equations 
(11151) and ( 11171) . In about three quarters of space, the root of the original 
exhausting function is negative, viz. for all negative values of z p and, if y p is 
sufficiently large, for positive values of z p too. The source lights the major 
part of space. 

By contrast, the mirrored exhausting function must be defined as 



A 



sm 



\J{x p ~ x q y + (y/(y p - y.y + zj + ^(y q - y_)* + z\f 



x p - x q ) 2 + (y p - y q ) 2 + (z p + z q f 



A | -|V A sm | if y P <y- + (y_ - y q )z p /z q with z p < 0, , 12Q , 
\ +| VA sm | elsewhere. 

Almost everything follows from (I119p replacing z q with — z q . The only excep- 
tion is the assignment of signs of the root. Using the freedom mentioned in 
the discussion behind equation fl65l) . your author reversed it. In about three 
quarters of space, the root of the mirrored exhausting function is positive, 
viz. for all positive values of z p and, if y p is sufficiently large, for negative 
values of z p too. Shadow prevails. 

Yet immediately over and under the screen, the signs of original and 
mirrored functions are the same, lim A sq = lim A sm , if the probe approaches 
the screen, z p — > 0; it does not matter if the approach takes place in inner or 
outer space, z p > or z p < 0. 

With these functions it is straightforward to construct solutions that sat- 
isfy the boundary conditions fl36l) exactly. Instead of (11051) we obtain 



. _ exp(ikR pg ) erfcy-z/cAsg exp(ikR pm ) ericy/-ikA sm | n/ ,_^ 
ak{r p ) — - - - - h U{K J 

it pq Z -LLpm Z, 

(121) 

exp(ikR pq ) erfc^-tkA sq exp(ikR pm ) eTfcy/-ikA sm 
h(r p ) = y — + + 0(k ). 

it pq Z flpm Z, 

(122) 

The symbol 0(k x ) was retained to remember that the differential equations 
(|41|) are not exactly fulfilled. 

The second terms on the right-hand sides of equations (11211) and ()122p 
considered isolated appear fantastic. They describe ghostly radiation from 
a source at r m that becomes bright only when it permeates the screen. Yet 
when they are considered in cooperation with the first terms, it is recognized 
that they describe the unavoidable reflection that is diffracted in a similar 
way as the primary wave. The solutions (I12ip and (11221) are valid in entire 
space. 
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From these solutions it follows that radiation emitted at r q and diffracted 
by a screen causes a second singularity at r m . In other words, diffraction 
creates a sharp image. Of course, the image is as weak as the multiplicative 
error function indicates, but it should be observable because the singularity 
sticks out. 

Folks might feel this prediction as daring. They might not appreciate the 
immense power of the methods developed here. The approximate solution 
(II 2 111122]) contains, as a limiting case, the only exact solution of diffraction 
problems known so far, namely Sommerfeld's celebrated stringent solution 

Sommerfeld's solution deals with the diffraction of a plane wave. Using 
the definitions in ( 15"5"|) and (11181) . R pq = \r p — r q \ and R pm = \r p — r m \, the 
spherical waves in front of the error functions in ( 11211) and (I122p can be 
expanded as 

g*Pg%^ = ^ fc | r «l) exp U k Zl±r p ) + 0(\r q n (123) 
I p q\ n?l ^ \^q\ ' 

exp(ifc|r p -r m |) = expggrj) ^ / \ + 2) 



The first factors on the right-hand sides do not depend on the position of 
the probe r p . They are absorbed in the normalizations of (I12ip and (I122p . 
When the source is moved to infinity, |r m | = |r g | — > oo, only the plane waves 
remain. Their fronts are defined by the normal vectors 

c q = lim - — f c m = lim - — y (125) 

|r q |->oo |r ? | |r,|^oo |r ? | 

with directional cosines as components 

C q G X C X ~t~ ^yCy ~t~ ^z^z C m ^x^x ~t~ ^yCy ^z^z • (126) 

Likewise the exhausting function for the plane-wave problem is defined 
as the limit 5 s (y) = lim A s (y) for |r c/ | — > oo. The straightforward calculation 
transforms (11151) to 



$s(y) = \jc 2 y + c 2 zS J (y p - y) 2 + z 2 - c y (y p - y) - c z z p . (127) 
Taking the value on the edge y = ?/_ yields the analog of (11191) 



S sq = ^Jc 2 y + c 2 ^ (y p - yJ) 2 + z 2 - c y (y p - yJ) - c z z p 

if y p <y^ + z p c y /c z with z p > 0, 
elsewhere. 
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:i28) 



and the mirrored exhausting function as analog of (I120p 



L = ^c 2 y + c 2 z^{y P -y-) 2 + 4 ~ °y(yp -v-) + c * z v 




-\\F^n\ if y p <y-- z p c y /c z with z p <0, 
+ 1 Vfisml elsewhere. 



(129) 



Thus the analogs of (jT2Tj) and (1X221 are 




sm 



(130) 




erfcy/^z 



(131) 



These functions form an exact description of diffraction. Not only the bound- 
ary conditions are satisfied perfectly, also the Helmholtz equations (HT1) as 
direct checking shows. One may be surprised at this, but it becomes com- 
prehensible when it is noticed that in this problem with the infinite edge 
and the primary plane wave, the only length is provided by the position of 
the probe r p . One can introduce new variables of position instead of r p by 
scaling the latter with with the inverse of the wave number k. Hence if the 
solution is known for one wave number, here for k — > oo, it is known for all 
wave numbers. 

Based on this argument, one supplements the proof of the fact that the 
two integrals in f lo^]) and f l5^|) are equal to T-2-7T / z p R pq , respectively. There 
is no edge at all. The only length is the distance between source and probe. 
One evaluates the integrals for k — > oo using stationary phase according to 
section |8] and introduces thereafter scaled variables to show that the result 
does not depend on k. 

Sommerfeld's solution is a special case of (ll30H13ip for flat incidence, 



i.e. c x = thus Jc 2 + c 2 z = 1. In this case, bkij p ) is proportional to the 

magnetic field Bfc(r p ) = —e x k 2 bk(r p ). The factor — k 2 is swallowed by nor- 
malization. In the problem with the other polarization, the electric field 
Efc(rp) is proportional to e a ,a/ c (r p ). Therefore Sommerfeld got along with- 
out the representation theorem of section [2J Moreover, Sommerfeld neither 
knew the method of stationary phase for two-dimensional integrals nor did 
he use the standardized complementary error function although it does not 
take much work to show that Sommerfeld's function is equivalent. Sommer- 
feld found his function via an ingenious contour integration, an approach 
which does not seem to admit generalization. Finally, it played a role that 
Sommerfeld was Felix Klein's pupil. Klein was the most influential advo- 
cate of uniformization meaning that roots had to be avoided because of their 
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ambiguities and to be replaced with suitable parametrizations necessitating 
transcendental functions. So Sommerfeld introduced, instead of the roots in 
(11281) and fll29j> . artifical angles, sines and cosines. Probably this rootopho- 
bia is the reason why Sommerfelds exceedingly important solution - it is the 
'harmonic oscillator' of diffraction - was never sufficiently appreciated and is 
not presented in most modern textbooks. 

Though Sommerfeld's solution holds only for flat incidence c x = 0, it is, 
due to translational invariance, easily generalized for skew incidence. This 
was done long before this article was drafted [31 § 11.6]. Moreover, the solu- 
tions f)130p and (I13ip can be conceived as Fourier components and used to 
build the most general diffraction at the straight edge. For example, diffrac- 
tion of a spherical wave, which is described by equations (I12ip and (I122p only 
asymptotically, was calculated by Macdonald [33]. Yet Macdonald's formu- 
lae are so difficult to survey that imaging by diffraction might have escaped 
notice. The present approach is simpler and generalizable. 

12 Diffraction by a slit 

Let us clarify next one of most lugubrious offences of traditional theory: the 
difference between Fresnel and Fraunhofer diffraction. To anticipate the an- 
swer, Fraunhofer diffraction is a misnamer. In the proper sense of diffraction, 
Fraunhofer diffraction does not exist. There is at best interference. It occurs 
only when a primary wave strikes an aperture with at least two opposite 
edges. The thus excited secondary waves interfere to generate a pattern now 
named Fraunhofer diffraction. The formulae derived in this article are pow- 
erful enough to describe the gradual transition from Fresnel diffraction to 
Fraunhofer interference. All is included in the equations (I103p and (I104p . 

To understand how the transition proceeds, the simplest example suffices: 
(i) Damping negligible; imaginary part — > 0; the wave number k is a 
positive number, (ii) Diffraction by a straight slit with edges at y± = ±yo', 
2yo is a positive number meaning the width of the slit, (iii) Normal incidence 
of a plane wave; thus c x = c y = 0, c z = 1, see (1125111271) : in equation (11041) . the 
plane wave replaces the spherical wave; consequently the exhausting function 
A s gives way to its descendant 5 S (I127p 



h(r p ) = exp(ikz p ) 



erfc J -ik5 s (y-) - eiic J -ikS s (y + ) 



+ 0(k~ 1 ) (132) 



2 



with 
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+ 



$s(y±)\ iiy p <±y , 
fis(y±)\ elsewhere . 



(133) 



(iv) The representative ak(r p ) is disregarded, (v) Effects of the reflected wave 
are neglected; z p > is implied. 

Equations fl 1 3 2 j) and fll33|> are still powerfull as the probe can be moved 
freely within the inner space z v > 0. Especially the equations remain valid 
when the probe approaches the screen z p — >■ 0. In most classical experiments, 
however, the probe is far away from the screen 
exhausting function 



z p — > oo. Therefore the 
can be expanded in terms of inverse powers of z p 



$s(yd 



y 2 P =F 2y p y + yj / | Q ^ y P Tyo ^ 



(134) 



It is worthwhile to mention that using the approximated exhausting function 
( 196]) yields the same result. The linearization (I94H95P spoils the approach 
to the screen. In the error functions of (11321) . which oscillate quickly, the 
first two terms on the right-hand side of fjl34j) must be retained because only 
the second term will yield specific results. But in algebraic expressions, the 
leading term will be enough, for example 



1 



8s(±Vo 



^(l + ft + 0((^ 



Vp 



(135) 



To pacify simple-minded scientists' dismay at the error functions in (11321) . 
let us return to elementary functions using the asymptotic expansion (1741) . 
Because of fl 133j) we must do it differently for y p — > — oo and y p — > +oo. For 
the latter case, equation ( 173]) has to be considered before ( 174"]) can be applied 



erfcJ-ik5 s (±y 



2 - 



exp(ikS(±y )) 



V-ink \J5{±y 



for y p +oo . 



(136) 



The minus signs from the last line of f!133j) and of the argument in ( I75p 
cancel. Opticians interprete this formula as representing the primary wave 
by the 2 and a secondary wave radiated from the edge by the exponential 
function. Therefore it is a formula for the domain of light. This seems to 
be a contradiction to physics since the domain of light is limited to \y p \ < y 
while we suppose here y p — > +oo. The contradiction is quashed since the two 
leading 2's cancel in 



Calculation based on (I134H136P yields 



—i exp ik(z p + 



Vv 



Zkyl sin(ky y p /z p ) 
2z p ')\\ Ti Zp ky y p /z p 



(137) 
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Notice the sine is brought about by the interference of the two terms in 
The formula is symmetric with respect to y p . It holds both for positive and 
negative values. A simple calculation for y p — > — oo along the lines just 
sketched confirms this. 

For a comparison with experiments, we need the time-averaged Pointing 
vector (II 10p . The probe usually is an absorbing layer spanned perpendic- 
ularly to the energy flux of the primary wave. The representative in (I137p 
produces via fl59l) and fl60l) with t = e x the electromagnetic fields 



B k (r p ) = -e x k 2 b k (r p ) E fc (r p ) = e y ukb k (r p )(l + 0(z p x )) + e 2 ... . (138) 

The dots behind e 2 substitute a function we presently do not need to know. 
The measured component of the energy flux is thus 

e 2 S(r p )~^|6 fc (r p )| 2 . (139) 

So = e z u)k 3 /(2/j,) is, because of the normalization chosen in (I132p . the energy 
flux of the primary wave. Straightforward evaluation of (I139p with (I137p 
gives 

s f \ is i o k y° { sin ( k yoVp/ z p)\ 2 f. An \ 
e z S(r p ) ~ | S 1 2^/0 — : T~^\ ' ( 140 ) 

irzp \ ky y p / z p J 

This is the most fabulous formula of traditional diffraction theory - Kirch- 
hoff's formula of Fraunhofer diffraction. However, Kirchhoff found only the 
factor with the squared sine. Historically, various normalizing factors were 
adjusted a posteriori, but the fabulous formula was never perceived as an 
equation among vectors. 

Here, by contrast, everything arises from first principles. The last two 
factors in f)140p when integrated over y p from — oo to +oo yield 1. The first 
two factors describe the primary energy flux squeezed through the slit of 
width 2y Q . Equation (11401) expresses energy conservation as warranted by 
Maxwell's equations in a medium without damping. 

The fabulous formula is worse than admitted in textbooks. From the 
estimate in (I135p we see it holds only for 



y « \y p \ « Zp . (141) 

yo << z p means that the probe is far away from the screen. This is all right 
because z p — > oo was announced at the start of the calculation. But the 
extendend condition (I14ip entails that Kirchhoff 's formula provides neither 
a valid description of the central peak around \y p \ < yo nor of large- angle 
diffraction \y p \ > z p . Only some intermediate wiggles are correctly seized. 
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That the condition y << \y p \ matters, is corroborated through the use of the 
asymptotic expansion (11361) . It is utterly wrong at values of the exhausting 
function close to zero. The triangle function and thus the exhausting function 
of optics is, according to (1641) . zero on the border of shadow defined here by 

2/0 = \y P \- 

The positive outcome of the preceding discussion is better understanding 
of the working of formula (I132p and its prototypes f 1 1 U 3 j) and (11041) . Let us 
position the probe first close to an edge. To describe this measurement, one 
of the terms in f 11 3 2 j) is enough. We must, however, utilize the error function 
full-fledged. When the probe is removed from the edge, but kept close to the 
screen, we are still in the realm of pure diffraction. Still one term of (1 1 3 2 1) 
suffices, but we may use now its asymptotic expansion (174"]) . Only in that part 
of space where distances to the edges are about equal, interferences show up. 
Both terms in (I132p must be kept. When the distances become large, we may 
use asymptotic expansions, but not at the borders of shadow. With increasing 
distances from the edges, the diffracted beams spread. In particular this is 
true close to the borders. The central beam is enclosed between these borders. 
Therefore it is always questionable to decribe diffraction and interference in 
the vicinity of the central beam using the asymptotic expansion f[74"|) . 

13 Diffraction by a circular aperture 

After it is confirmed that the new theory encloses correct parts of the tradi- 
tional theory as special cases, we are free to step into the entirely unknown. 

The round stop is the most often installed component in optical instru- 
ments. We want to calculate the diffraction it causes, more precisely its 
Fresnel diffraction, since some information on Fraunhofer interference is al- 
ready known - Airy's formula [31 § 8.5.2]. The edge of a circular aperture is 
not straight. Nevertheless the universal formula of diffraction (11051) applies 
again. All we have to do is to find the suited exhausting function. 

The suitable coordinate system is the cylindrical one introduced in section 
El see dSB. 

- 7r < if < +7T 0<p<p = const z = (142) 

defines a circular aperture in the screen. The exhausting dependence is found 
from equations (1551) and (1561) 

d v R p {p, f) = -d v R q {p, f) (143) 

with 

Rp(p, f) = ^]p 2 p + p 2 - 2p p pcos(if p -f) + z 2 p 
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RM <p) = \jp\ + p 2 - 2 p,p C0S K -<p) + z l • ( 144 ) 

For normal incidence, p q = 0, one can observe the solutions using the rubber- 
ribbon experiment explained in section [SJ Actually, there are two: 

ip = (fp and if = ip p + 7r . (145) 



These exhausting dependencies <p = analogues of £ = H s (t7) in 

surprise as they do not depend on p. Inserting them into the triangle function 
(l56|) gives the exhausting function of the circular stop. We need it at the edge 
P = Po- 



A, = \f(Pp ~ Po) 2 + zl + Jfi + z\ - y/fi + (z p - z q f . (146) 
The definition has to be completed with the disposal of the root 



/A7| if p p > p {z q - Zp)/z q with z p > 0, 
/AJ elsewhere, 



(147) 



according to equations (165 j) and (|68p . The exhausting function of the second 
solution in f)145p is derived from this replacing p with — p . Your author 
names the first solution with the positive po near because its point on the 
edge is closer to the probe than in the opposite case. The second solution 
with the negative p is thus called far. 

For mathematical reasons, only the near solution may dominate. The 
definition ( I142p of cylindrical coordinates restricts the range of the angle ip. 
Yet unrestricted variation of the variable £, analog of (p, is a premise of the 
principle of utter exhaust as declared in section |HJ The violation of this 
premise is not disastrous as long as the exhausting function A(£, 77) takes a 
sharp extremum on the edge. In this case, the finiteness of integration causes 
negligible corrections. But for the round stop A(£,r?) becomes equal for all 
points on the edge when source and probe both approach the optical axis; 
there is no extremum of A at all. Consequently the universal formula 
of diffraction f ll05p applied to circular apertures needs p p > z p > to secure 
its applicability. 

Of course, the primary wave is also diffracted at the far side of the circu- 
lar stop. There will be similar interferences as those caused by the slit, see 
section [121 but the contributions from the far side are here, in the range of va- 
lidity just declared, small. Fraunhofer interferences cannot be described well. 
Nevertheless, within the restricted range, there seems to be no rival of the 
formulas f)146p and f ll47p . We can calculate diffraction immediately behind 
a circular aperture and predict what a curved edge does to polarization. 
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For general parameters, the exhausting equation (11431) cannot be solved 
in terms of roots. When the transcendental functions sine and cosine are 
rationalized introducing, for instance, the variables c and s 

c = cos(</? ) s = sm[ip ) (148) 

where c 2 + s 2 = 1, and when the roots are removed by squaring, an algebraic 
equation of the sixth degree is obtained. 

Nevertheless we can find more physically relevant solutions, even for 
oblique incidence, performing the rubber-ribbon experiment as above. Namely 
the one-sided Hat situation where source and probe sit on the same side of 
the optical axis, and the two-sided Hat situation. In both situations, there 
are near and far solutions: 

one-sided flat: p> p = p> q near: ip = ip p far: ip = ip p + ir (149) 
two-sided flat: <p p = p> q + n near: ip = (p p far: (p = (p p + n . (150) 

But this is not enough when vector fields as in fll08[) and (11091) are to be 
calculated. For them we must differentiate the exhausting function with 
respect to z p , p p and (p p . The last differentiation is not feasible when the 
function is only known for a fixed (p p . 

Mending the shortcoming is easy. The sine S and the cosine C were 
defined in (|69i) . In the context of (I143p . they must be considered as given 
and fixed. The sine s and the cosine c were defined in (11481) . They contain 
the angle which is sought as the exhausting dependence <p = $ s (p). With 
the abbreviations 

P( P ) = ^{p p - pf + zl/ Pp Q(p) = ^{p q - pf + zl/p q (151) 

the exhausting equation (I143p can be written in four equivalent variants 

Cs-Sc Cs + Sc 

yJP{±py-2p{Cc + SsTl)/p P ~ ^Q(± P ) 2 -2p(Cc-Ss^l)/p q ' 

(152) 

Differently to all other equations in this article, signs can be chosen indepen- 
dently. Just under the roots, the signs must be altered coherently. 

These variants permit to find solutions in the neighborhood of (I149p and 
(11501) such that they are exact in linear terms. From (|69|) we have 
one-sided: S = 0(p> p — <p q ) C — 1 + 0((p> p — <p q ) 2 ) 

two-sided: S = 1 + 0((p> p — <p g — ix) 2 ) C = 0(ip p — <p q — 7r) . 
Inserting the estimates according to (1149111501) 
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one-sided near 
one-sided far 

two-sided near 
two-sided far 



ip = <£p + 0((p p - (p q ) 

(f = (fp + IT + 0(lf p ~ f q ) 

(f = ip p + 0(<p p -ip q -n) 

V? = ip p + 1T + 0{if p -ifq- 



7T 



in the definitions (I148P yields 



one-sided near 
one-sided far 

two-sided near 
two-sided far 



0(tp P - ip q ) 
0(<p P - <Pq) 
l + 0(( Vp - 



Vq ~ 71 ) 



1 + 0((<p p - <p q )*) 
-1 + 0(fo P - fft) 2 ) 

0(cp p -(fq-TT) 
0(ip p -(fg-w) . 



S 

s 
s 

S = -l + 0{(ip p -(fq- 7r) 2 ) 

Hence the nominators in (11521) are always small of first order. In the de- 
nominators, -P(±p) and Q{±p) are positive constants. But the expressions 
Cc+Ss =p 1 and Cc — Ss =F 1 are always small of second order if signs are suit- 
ably selected. Thus these expressions can be neglected. Without extertion 
one finds now the exhausting dependencies for the one-sided flat situation 



near: 



-S 



p(p) - Q(p) 



p(p) + Q(p) 

and for the two-sided flat situation 

P(p) + Q(- P ) 



far: 



S 



P{- P )-Q{-P) 
P(-p) + Q(-p) 



near: 



-c 



p(p) - Q(-p) 



far: c = C 



P{-p) + Q{p) 
P{-p)-Q{p) 



(153) 



(154) 



The exhausting function of the one-sided flat situation, based on the near 
solution and taken at the edge p = po, is, when abbreviations are restored, 



A s = J (pp - Po) 2 + 4 + J(p q - po) 2 + z\ 



~y(Pp- Pq) 2 + ( Z P~ Z q 



+2 



(.Pp-POY + Zl J (Pq - PO) 2 + Zq , 



PpPO PqPO 
V '(PP ~ Pq) 2 + ( Z P ~ Z P) 2 s-l\ . 2^P~Vq 



PpPq 



sin 



(155) 



/ s:=+Ia/s:i if 



pp > 



Po(Zq - Zp) + p q Zp 



1 + 2- 



Pq Z p . 2 ( Pp~ i Pq* 



sm 



with z v > 0, (156) 



Z q Po\Zp — Z q ) I 

\f~A s = — | v^II elsewhere. 

The equations (I68H69P were used to determine the border of shadow. 
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The respective expressions the two-sided flat situation are 



A* = J(Pp - Po) 2 + *% + J(pq + Po) 2 + 



-\j{pp + Pq) 2 +{Z P ~ Z q f 

t \/(Pv~ f. 

+2 



(p p -po) 2 + Z 2 p y \ Pq + po) 2 + Z* 



PpPO PgPO 



. v (Pp + Pg) + (z P ~ z p) >_A 2 y P -y g , * 

+(- cos g 157 

/SJ=+|VSJ| if 

Pp > Po( " 9 ~" p) ~ P ^ (l-2 cos 2 with z p > 0, (158) 

Zq P0\Zp — Z q ) 2 

/AT = — l-v/ATl elsewhere. 



All exhausting functions based on the far solutions can be obtained from 
these formulas inverting the sign of po- The errors are 0((ip p — <p?) 4 ) or 
0(((fi p — ip q — vr) 4 ), respectively. 

Generally, one may solve equation (I152p numerically. The effort for this 
is orders of magnitudes less than for a direct numerical solution of Maxwell's 
equations. 



14 Outlook 

As we all know, everthing was foretold by Aristotle or Jesus Christ. So it is 
only with timidly throbbing heart that your authors dares to ask the deciding 
question: What might be novel in this article? 

All monographies and textbooks pretending to teach the solution of Max- 
well's equations should be committed to paper recycling. The same applies 
to those treatises wherein optics is allegedly derived from electrodynamics. 
The general representation theorem in section [2] beats them all. That the 
boundary conditions on metallic surfaces simple for the representa- 

tives as deduced in section El this was never anticipated. These discoveries 
are synthesized with Sommerfeld's criticisms and proposals in the sections H] 
and [5j The thus derived solution of electromagnetic diffraction polarization 
inclusive, described in section O is a new finding. Compared to this, the 
discussion of the triangle function in section O appears at first sight of minor 
importance. But the root of the triangle function and a consistent assignment 
of its signs are of prime importance for applications. The big thing, however, 
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is the principle of utter exhaust declared in section [SJ For the first time, 
stationary phase is applied where it is really needed. Utter exhaust seems to 
have no predecessors at all. It is so powerful, it would be explained in dozens 
of textbooks if it was known. Utter exhaust shows its power in the universal 
formula of diffraction, see section |9j Its application to diffraction by an edge 
in section [TD] yields a comprehensive formula which encloses, as a limiting 
case, Sommerfeld's stringent solution, see section [TTJ The great wonder is 
the perfect agreement though utter exhaust is just an asymptotic evaluation 
for short wavelengths. For diffraction by a slit, as discussed in section Q21 
utter exhaust copes with Fresnel diffraction and Fraunhofer interference in 
a unified way. The gradual transition between these fundamentally different 
patterns can be calculated. The derivation of these formulae and their no- 
tation is new. Everything is now easier to survey than in previous theories 
The formulas in section [13] describing diffraction by a round stop, derived via 
utter exhaust, are completely new. 

Hundreds of physicists will find rewarding foundations for future busi- 
ness. The physical discussion of the formulas derived here has only begun. 
Important systems, e.g. diffraction by a hook or a cusp, were not yet an- 
alyzed though they are now accessible. Moreover, the circular aperture is 
an especially rich system which deserves further studies. The foundation of 
its theory is the integral f l58|) . Here it was evaluated for short wavelengths 
only, but more asymptotics must be considered. Evaluation for large dis- 
tances from the screen will yield Fourier integrals and thus comfort all those 
who understand imaging as a sequence of Fourier and inverse Fourier trans- 
forms. The asymptotics for small angles are scarcely understood though they 
are most important for applications in optical systems. Quantitative com- 
parisons with experiments are altogether rare. The best measurements are 
probably those made with microwaves [31]. 

At last we have a reasonable theory of electrodynamic diffraction, but 
also part of acoustics will change. The methods developed from section [3] on 
apply to the velocity field of sound if the representatives dk(r p ) and bk(r p ) are 
construed as scalar potentials. The boundary conditions derived in section 
[3J are just the most often used ones in acoustics, namely on the soft or the 
hard wall, i.e. acoustic impedance zero or infinite. Almost everything remains 
the same, in particular the universal formula of diffraction f ll05p . Only the 
dispersion relation of sound is different from (T421) since compressible fluids 
refuse the telegraph equation [17] . 

It is time to reconcile with Gustav's fans. Your author knows the pit- 
falls of pioneering. He reveres Gustav Kirchhoff affectionately and repents 
bitterly not to have read his books first. Almost all epigone scientists copied 
Kirchhoff 's ideas, and it is sometimes felt that the copies were made without 
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understanding. In other words, modern scientists do not attain KirchhofT 's 
intellectual level. It is, however, remarkable that just Max Born copied Kirch- 
hoff's. Born is reputed as the founder of quantum mechanics in so far as he 
invented the probabilistic interpretation. A careful comparison of Born's 
monography |2] with Kirchhoff's textbook on optics [6] revealed that Born 
copied the probabilistic interpretation, too. Kirchhoff, an adherent of Ger- 
man idealistic philosophy according to which material objects descend from 
ideas, made a theory for Verriickungen, i.e. for nonobservable items. The 
reader might remark the contrast to Maxwell who based his theory on mea- 
surable quantities, namely on the electromagnetic fields. Nevertheless Kirch- 
hoff, being also a physicist, desired to compare his ideas to measurements. 
So he stipulated, without any reason, that the squares of those Verriickungen 
were intensities, i.e. nonnormalized probabilities: "Von dieser hangt die In- 
tensity J des Lichtes ab. Als Maass derselben fiihren wir den Mittelwerth 
des Quadrates der Verruckung eines Aetherteilches ein" [61 p. 9]. Kirchhoff 
is the true inventor of the probabilistic interpretation, but the invention was 
not a sign of mental strength, but lack of ability to solve the correct differ- 
ential equations correctly. As shown in this article, in particular in section 
[T2| electrodynamics do not need a probabilistic interpretation. 

The Dirac equation, actually a system of coupled partial differential equa- 
tions, does the same for electronic motion as Maxwell's equations do for pho- 
tonic propagation. A probabilistic interpretation is not necessary because the 
Dirac equation describes observables, namely the external electrical charge 
and current densities which were discussed in section [2j The Dirac equa- 
tion even defies that interpretation. Eighty years after Dirac's discovery, no 
compatible probability formalism for the Dirac equation is known. 

Your author meanwhile derived representation theorems for the Dirac 
equation [32] • They simplify its solution in the same way as the represen- 
tation theorem of section [2] simplifies the solution of Maxwell's equations. 
Schrodinger's ip turns out to be a representative similar to the functions a 
and b used here. The Schrodinger equation appears as an auxiliary simi- 
lar the telegraph equation. Interpreting Schrodinger's mathematical tools 
physically appears as Verruckung again. 

Moreover it is already known that several effects seeming to prove prob- 
abilistic behavior in quantum electrodynamics can be derived instead from 
the spatial distributions of charge and current as determined by the Dirac 
equation [37]. Nevertheless, a certain group of experiments exhibits random 
behavior of photons and electrons, viz. absorption by small absorbers [38J. 
Yet the respective dynamical problems were never solved - as pipe turbulence 
was not solved before 1988 [27J . Random behavior in quantum electrodynam- 
ics might stem from deterministic chaos similarly as in other areas of nature. 



43 



It might happen that Heisenberg's mysticism regarding quantum "particles" 
and "uncertainty" will be gone soon. 
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